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Dust aerosol plays a fundamental role in the behavior and evolution of the martian atmosphere. The first 
five Mars years of Mars Exploration Rover data provide an unprecedented record of the dust load at two 
sites. This record is useful for characterization of the atmosphere at the sites and as ground truth for orbi- 
tal observations. Atmospheric extinction optical depths have been derived from solar images after cali- 
bration and correction for time-varying dust that has accumulated on the camera windows. The record 
includes local, regional, and globally extensive dust storms. Comparison with contemporaneous thermal 
infrared data suggests significant variation in the size of the dust aerosols, with a 1 pm effective radius 
during northern summer and a 2 pm effective radius at the onset of a dust lifting event. The solar longi- 
tude (L s ) 20-136“ period is also characterized by the presence of cirriform clouds at the Opportunity site, 
especially near L s = 50° and 1 1 5°. In addition to water ice clouds, a water ice haze may also be present, 
and carbon dioxide clouds may be present early in the season. Variations in dust opacity are important 
to the energy balance of each site, and work with seasonal variations in insolation to control dust devil 
frequency at the Spirit site. 

© 2014 Elsevier Inc. All rights reserved. 


1. Introduction 

Dust aerosol plays a fundamental role in the behavior and evo- 
lution of the martian atmosphere. Atmospheric dust and ice clouds 
have a direct effect on both surface and atmospheric heating rates, 
which are basic drivers of atmospheric dynamics. In turn, the 
resulting atmospheric motions influence the distribution of the 
dust itself, thus setting up a complex feedback {e.g., Newman 
et al., 2002). For example, the lifting of dust from the surface of 
Mars is a function of winds and atmospheric instability (i.e., turbu- 
lence) at and near the surface, both of which are influenced in turn 
by heating rates dependent of the amount of atmospheric dust 
loading. This type of non-linear feedback mechanism is suspected 
of determining whether in a given year Mars will see a global dust 
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event or only an assortment of local to regional storms (e.g., Rafkin, 
2009). In addition, the changes in the dust column can directly 
influence the vertical scale of the atmosphere, and sometimes pro- 
duce large variations in density (observed in the form of aerody- 
namical drag) at orbital and aerobraking altitudes (Bell et al., 
2007). 

The importance of dust was recognized early in the exploration 
of Mars (Gierasch and Goody, 1 972) and consequently dust became 
a target of atmospheric studies (see review by Smith, 2008). While 
earth-based and orbital observations comprise a large part of the 
historical record, landed observations of optical depth represent 
a critically important component because of the higher accuracy 
and precision generally available from such measurements. That 
is to say, although orbital observations allow for a more complete 
characterization of dust loading over spatial and temporal scales, 
they generally embody more modeling or retrieval assumptions 
than one might find with surface observations (e.g., Smith, 2004, 
2009; Wolff et al., 2009). In essence, they cannot provide generally 
the self-consistent “ground-truth”, the proverbial “grail” of remote 
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sensing algorithm validation. Prior to the Mars Exploration Rover 
(MER) mission, the record of surface observations comprised about 
1.3 martian years of data at each of the two Viking Lander sites 
(Colburn et al., 1989) and about 0.1 martian year of data at the 
Pathfinder site (Smith and Lemmon, 1999). However, despite this 
temporal coverage, there were no simultaneous orbiter/ground 
station observations. In fact, even the application of these data to 
close comparisons has been fairly limited (Hunt, 1979; Martin, 
1986; Clancy and Lee, 1991). 

In January 2004, a new era of surface-based atmospheric obser- 
vations began with the arrival of the MER missions. The two rovers, 
MER-A (Spirit) and MER-B (Opportunity) landed shortly before the 
Southern autumn equinox on opposite sides of the planet in Gusev 
crater (14.57°S, 175.48°E, -1.9 km) and in Meridiani Planum 
(1.95°S, 354.47°E, -1.8 km), respectively (Arvidson et al., 2006; 
Squyres et al., 2006). Although “mission success" was considered 
to be 90 sols (a sol is a martian solar day, about 1.0275 Earth solar 
days) of operation for each rover, Spirit functioned through sol 
2209 and Opportunity has accumulated more than 3500 sols as 
of this writing. Throughout their mission, the rovers have provided 
atmospheric optical depth measurements though direct solar 
imaging (Lemmon et al., 2004). The evolving optical depth data 
set has been used in coordinated (“ground truth") and contextual 
observations of the atmosphere and surface (e.g., Wolff et al., 
2006, 2009, 2010; Mateshvili et al., 2007; Clancy et al„ 2007; 
Vincendon et al., 2007, 2009; Heavens et al., 2011; Hoekzema 
et al., 2011; Petrova et al., 2012; Montabone et al., this issue). In 
addition, this dataset has been applied to investigations of surface 
activity (e.g., Arvidson et al., 2006; Greeley et al., 2006a, b; Rice 
et al., 2011), photometric models investigating surface or atmo- 
spheric properties (e.g., Johnson et al., 2006a, 2006b; Bell et al., 
2006b; Kinch et al., 2007; Shaw et al., 2013), and as inputs to dis- 
cussions of chemical properties or biopotential (e.g., Schuerger 
et al., 2006, 2008; Vincendon et al., 2010). Although the MER solar 
observations were briefly described by Lemmon et al. (2004), the 
scope and the complexity of the current dataset have greatly ex- 
ceeded that presentation. In this paper, we describe in detail the 
generation of MER optical depth record and the associated uncer- 
tainties, as necessary to understand the archival data. Section 2 
describes the observations, their operational context, and the 
on-board processing done by the rovers. Section 3 describes the 
processing of the downlinked data into an optical depth archive. 
Section 4 presents a discussion of implications of the data set, 
and Section 5 presents the conclusions. 


2. Observations 

2.1. Overview and rover operations context 

The MER atmospheric optical depth archive spans more than 
2200 sols for Spirit and more than 3500 sols for Opportunity. The 
mission began near the end of Mars year (MY) 26, using the Clancy 
et al. (2000) convention for numbering. In this work, we consider 
the entire Spirit data set and the Opportunity data set through 
the northern spring equinox at the end of MY 31 (31 July 2013). 
Pancam solar images using two neutral-density (ND) coated filters 
form the fundamental measurements. We first describe the imag- 
ing and raw Experiment Data Records (EDRs) in this section; de- 
tails of the optical depth retrieval are provided in the next section. 

For context, the overall pattern of optical depth is shown in 
Fig. 1, and is not strongly sensitive to the calibration details de- 
scribed below for sols 1-1200. Seasonal context is discussed in Sec- 
tion 4, but generally low and stable opacities signify southern 
winter (and aphelion), while high and variable opacities correlate 
with southern summer (and perihelion). 


Pancam has two cameras (eyes) per rover, each with a CCD 
detector, a filter wheel, and optics that map a 16° x 16° field of 
view onto 1024 x 1024 pixels (Bell et al., 2003). The solar filters 
are filter number eight in each of the right and left Pancams. The 
nominal wavelengths are 440 nm (L8) and 880 nm (R8); compara- 
ble Viking observations were done at 670 nm (Colburn et al., 1989), 
and Pathfinder observations were made at 440, 670, 880, and 
990 nm (Smith and Lemmon, 1999). Late in pre-launch calibration, 
a red-leak— excess long wavelength sensitivity— was discovered in 
the L8 filters. This defect was not significant prior to application of 
the ND coating, and was imprecisely measured after the coating 
due to low signal. The leak affects a broad wavelength range, 
roughly centered at 750 nm. The red leak is somewhat more signif- 
icant in the pre-flight calibration, using a cooler lamp, than in 
flight, imaging a hot Sun. Nonetheless, in general we refer to the 
optical depths by filter name, rather than wavelength. 

Most solar images were taken as part of a dedicated campaign 
to create a detailed optical depth record. Additional images were 
also obtained in order to characterize the optical depth as context 
for another primary observation, such as sky or surface photome- 
try. In addition, the Sun’s apparent motion over time was used 
by software aboard the rovers to determine spacecraft attitude 
after some drives. A subset of the latter (“engineering”) images is 
included in this analysis. A sequence name of P26xx (e.g., P2601) 
indicates an atmospheric-science sequence; sequences named 
P27xx was used for some lunar transit imaging; and F0006 desig- 
nates an “engineering” sequence. An example product identifier is 
2P127616614ESF0309P2600R8M1, where ’2P' indicates Spirit's 
Pancam, there is a 9-digit time stamp, a product type (ESF = sub- 
frame, ETH = thumbnail image), a 4-character location identifier, 
sequence (P2600), filter (R8), a producer code and version number. 
The product identifier is included in the archival record described 
in the next section. 

One of the standard monitoring sequences was routinely done 
shortly before noon (specifically immediately before or after the 
transition from one sol’s master sequence to the next) to maximize 
completeness of a record for monitoring performance of the solar 
arrays. All solar images were exposed for a commanded duration 
rather than using auto-exposure; imaging when the Sun was low 
or the optical depth was high required longer exposures that are 
typical of some (e.g., P2633) sequences. Observations were 
commanded generally as paired L8-R8 exposures. Exceptions in- 
clude some Phobos and Deimos transit observations (R8 only) 
and F0006 sequences (L8 only). A somewhat larger subset of the 
observations in the record are not paired. Reasons include: differ- 
ences in exposure times resulting in one image saturating or 
under-exposing when conditions rapidly change; failure in the on 
board Sun detection, sometimes due to a faint Sun and a bright cos- 
mic ray; and unrecovered loss of data during transmission. 

Fig. 2 shows the time-of-sol sampling for each rover. The nom- 
inal operational day for the rover was initially about 1100 to 1700 
HLST (hybrid local solar time— this is the local mean solar time, off- 
set to make it close to local true solar time during the nominal mis- 
sion). On certain sols, imaging earlier and later was possible, 
particularly when solar power levels were high (e.g., Southern 
summer and/or clean solar arrays). Early-morning imaging was 
especially rare due to the low temperatures and the need to heat 
actuators and electronics prior to use. In addition, frequent optical 
depth observations were made early in the mission; as mission 
longevity raised concerns about actuator lifetime, the number of 
observations per sol was scaled back. Patterns seen in Fig. 2 reflect 
the operational part of the sol and the tendency to fix opacity mea- 
surements to communications windows: direct to/from Earth com- 
munications near 11 HLST, and Mars Odyssey overflights near 17 
LTST (local true solar time, defined by local true midnight, and 
varying in sol duration due to Mars’ eccentric orbit). The figure also 
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Fig. 1. The MER optical depth record for Opportunity and Spirit (offset by 2). Opacity observed in L8 (blue symbols) and R8 (red line) is shown. Mars year is 27.0 at northern 
spring equinox (L s = 0°, Opportunity sol 41 ) following landing and increases proportionally to L s between successive northern spring equinoxes; minor ticks indicate 45° of L s . 
A sol scale along the top axis indicates Opportunity sol number (20.5 less than Spirit sol number), which is non-linear in L s . All valid measurements are presented. (For 
interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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Fig. 2. Local true solar time (LTST) for each optical depth measurement for (A) Spirit and (B) Opportunity. LTST midnight is defined as 0 h and 24 h when the Sun is at its 
lowest ephemeris position and increases linearly with time through a sol. The pattern near 12 h is associated with the direct to/from Earth communication; patterns near 16 h 
are associated with Odyssey relay communications. 


illustrates the changing relation between HLST and LTST (the 
equation of time) in the pattern oscillating around 12 LTST. 

The procedures and data products described in this paper were 
developed in conjunction with evolving operations concepts. The 
amount of data collected and the specifics of the implementation 
varied with lessons learned and in response to operational pres- 
sures. A key focus of the process was to keep a quick-look optical 
depth available on a rapid timescale, while maintaining consis- 
tency with the final product. 

2.2. Rover on-board processing and downlink 

A typical solar-image command used the rover’s on board 
ephemeris to aim the camera at the Sun. This approach allowed 
for flexibility with respect to the time a sequence is run (e.g., 
P2600 was run at many times of sol, independently of rover 
attitude or the Sun’s position), but was susceptible to drift in the 
Sun’s position in the image frame when a rover’s attitude self- 
knowledge is degraded. Attitude knowledge was typically good 
to within 3°, and was restored with dedicated solar imaging 
(F0006 images) and Mini Inertial Measurement Unit (M1MU) use 


based on operational needs. The typical solar-imaging command 
also resulted in 4 EDRs being returned to Earth, a subframed image 
and a thumbnail image in each of the two filters. "Sunfind”, or 
auto-subframed, images, are 63 x 63 pixels from the original im- 
age, centered on the ~20-pixel diameter Sun; manual subframes 
were required to be larger parts of the original image, to account 
for pointing and timing uncertainties. Thumbnail images were 
automatically created for all MER images and “downlinked” with 
high priority for quick operational assessment of a sequence’s re- 
sults. These thumbnail images were available before the subframes 
(sometimes by weeks), and in some cases remain the only products 
available (e.g., due to failure of the automatic subframing or loss of 
data). 

Exposure durations were typically 3-5 s for L8 and 0.3-0.4 s for 
R8, but exposures of up to 5 min were possible. However, exposure 
times of over 1 min “trail” the Sun by at least its own diameter and 
were avoided. Pancam’s image acquisition relies on an electronic, 
rather than physical, shutter. That is to say, image readout is con- 
tinuous until the start of an acquisition, when readout pauses for 
the exposure duration and then a readout is recorded. A zero-sec- 
ond (“shutter”) exposure would see only a shutter smear caused by 
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the 5 ms readout while the detector is exposed to light. Longer 
exposures would superimpose the image on that smear. The smear 
could be corrected out at image acquisition by subtracting an 
actual zero-second image or after downlink by analytically synthe- 
sizing a shutter image. Due to the Sun’s apparent motion between 
actual exposures, the synthetic shutter method was typically relied 
on; operational issues prevented that mode from being used 
exclusively. 

To minimize data volume, solar images were subframed to be as 
small as practical. The automated subframing algorithm relied upon 
a digital filter applied to the stored image, looking for the maximum 
brightness of a central box slightly larger than the Sun after sub- 
tracting the average of the pixels adjacent to the box. If the summed 
data numbers (DN) exceeded some threshold, a 63 x 63 subframe 
was extracted and stored for downlink. In order to image the Sun 
when the predicted brightness might fail to meet the threshold, 
manual subframes were used. Such cases included low Sun and high 
optical depth situations when even long exposures would be faint, 
or when the expected variability was sufficiently high that the range 
of possible brightness was too large for the range between the 
threshold and detector saturation. Images were sometimes com- 
panded from 12- to 8-bits per pixel via on board look-up table and 
were generally losslessly compressed. 

Thumbnail images, 64 x 64 images that have been downsam- 
pled from 1024 x 1024 full frame images using pixel averaging, 
then compressed to 4 kbit each, were created for all MER images 
and downlinked at higher priority than any regular images. As 
such, the solar-image thumbnails formed the basis for a quick-look 
optical depth, being processed to optical depth in the same manner 
as the full-resolution images. In the cases where they are the only 
record, the thumbnail data are merged into the archive. 

A lesson for operation of future missions is that the metric for 
total DN within a local-background-subtracted sun-sized box is a 
meaningful number. This single number is roughly proportional 
to flux, subject to a few perturbations described below. The 
on-board processing was one step away from downlinking an 
operationally useful quick-look optical depth as a single number 
attached to an image label, with no actual image needed. The initial 
build of flight software returned this flux-like number for even 
failed sunfind attempts (e.g., when the Sun was detectable but 
failed to meet the threshold), thus allowing optical depth retrieval 
from the image header even when the conservatively set threshold 
caused rejection of an otherwise valid image. The intensity gap in 
which images could be processed to provide a useful optical depth 
measurement yet fail to meet the threshold is due to an operation- 
ally necessary compromise. Science use of the data would motivate 
a low threshold in order to gather the most data; however, 
engineering use required a high threshold in order to avoid false 
positives influencing attitude derivation, and only a single thresh- 
old was available. 


3. Optical depth calibration 

3.1. Overview 

Conversion of solar images into an optical depth measurement 
is fundamentally simple: following the Beer-Lambert law, 

F = F T oa ■ e~ T i, where F is the measured flux, F T0A is the top-of- 

atmosphere flux, z is the desired normal-incidence extinction opti- 
cal depth, and ;/ is the airmass, or ratio of slant-path optical depth 
to normal optical depth. Significant variations in the true F T oa var- 
ies over time inversely to the square of the known Mars-Sun dis- 
tance. Fig. 3 shows example solutions using subsets of the data. 
The trick is in the accurate determination of F and F TO a in self-con- 
sistent units to minimize systematic errors, and to a lesser extent 



Fig. 3. Beer-Lambert Law fit. Opportunity L8 fluxes, normalized to a Sun-Mars 
distance of 1.50 astronomical units, from sols 23 (filled circles), 179-184 (filled 
squares), and 1 691-1694 (open diamonds) are shown as a function of airmass along 
with exponential fits for each period. The early data has nearly the same 0-airmass 
intercept, but different opacity (0.90, then 0.47). The post-dusting data have a 
moderate atmospheric opacity (0.57) plus an additional 50% extinction from 
window-dust. Uncertainties are smaller than the symbols. 

in the determination of )/. The methods employed are described 
in this section. 

Derivation of atmospheric transmittance, FIF TO a, was done via a 
relative calibration approach similar to that in the Mars Pathfinder 
analysis of Smith and Lemmon (1999). Although the absolute cali- 
bration of the L8 and R8 filter response was sufficient to determine 
exposure times, comparison of an instrumentally calibrated flux 
with a ‘‘known’’ solar flux would result in errors that were large 
compared to the amplitude of interesting atmospheric variability, 
and specifically would introduce spurious diurnal variability con- 
trolled by airmass. Ideally, F toa would be directly observed, but this 
was impossible. Instead, one may take advantage of the linearity of 
ln(F) = ln( F 7 oa) — t ■ rj. By asserting that t is unchanging in a period 
containing several measurements at different the slope (t) and 
intercept (ln(F rcM )) can be retrieved. Complications associated with 
this methodology arise from systematic and random natural 
variability, and from the fact that (by the end of the observational 
record) the cameras had spent 9 years in a dusty environment. We 
discuss specific solutions below. 

Image calibration for optical depth retrieval was done via a ded- 
icated pipeline. Image calibration was done outside of the standard 
calibration pipeline (Bell et al„ 2006a) for the following reasons. ( 1 ) 
There was a mix of images with and without onboard shutter sub- 
traction, requiring those without to get synthetic shutter correc- 
tion for consistency. In general, subframed images cannot and 
did not have a synthetic shutter generated. In the special case 
when the subframe contains the brightest object around by orders 
of magnitude, the generality does not hold. (2) Calibration must 
subtract a dark current and bias value from the images. The stan- 
dard pipeline did not report negative intensities, effectively adding 
flux to the background when the subtraction was overdone. Accu- 
rate background correction is important, as described below. (3) It 
quickly became apparent during operations that control over the 
temperature and flat field corrections would help to produce more 
accurate results. 

3.2. Image processing 

A raw image experiment data record for a subframe was ex- 
panded from 8 to 12 bits via inverse lookup table when needed. 
The image was then calibrated using elements of the Pancam 
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calibration pipeline (Beil et al., 2006a), following the Bell et al. 
(2006a) flow from image input to flat-field removal. The procedure 
involved subtracting out bias and dark current and the tempera- 
tures and exposure time reported in the image header; subtracting 
a model shutter image if needed; then dividing out the exposure 
time. Thus, F y = (S(/j) - B(T) - D(i,j,T ) ■ t exp - Sh{ij))lt exp , where F u 
is the derived flux for pixel ij in DN/s; S is the observed signal in 
DN, B is the modeled bias at temperature T, D is the modeled dark 
current rate, and exposure duration t exp , and Sh is the modeled 0-s 
shutter or readout image. 

Flux extraction required finding the Sun as the brightest pixel in 
a boxcar smoothed version of the image (a small number of false 
positives were manually screened). The Sun’s position was then 
determined more precisely via the centroid in the original im- 
age-note that this position was almost always the image center 
for auto-subframed images. The background flux was determined 
in an annulus ~20-30 pixels (2-3 solar radii) from the center. 
The residual background came from a combination of error in bias 
and dark subtraction and the sometimes-detectable, forward-scat- 
tered diffuse light. The influence of scattered light increases with 
airmass, so the last step was needed to avoid systematic errors. 
After background subtraction, pixels were summed to retrieve a 
flux in DN/s and normalized to a constant Sun-Mars distance of 
1.5 AU. Images with saturated pixels were flagged with a negative 
flux and discarded. In addition, images containing transits of Mars’ 
moon across the Sun (Bell et al., 2005) were flagged and discarded. 
One step that we omitted from the standard calibration is the flat 
field correction, which is discussed in the next subsection. 

Determining flux from thumbnail images is conceptually similar, 
although only the image header is needed. On-board background 
subtraction removed contributions from dark current and detector 
bias, and integrated flux was returned in the header along with the 
location of the Sun in the original image. In the pipeline, saturation 
was flagged for a too-high integrated flux (determined by inspection 
of subframes with saturation), or a too-high maximum DN level in 
the thumbnail. Exposure time was then divided out. Such use of 
the header value is actually preferable to processing the downlinked 
thumbnail image due to the issue of lossy compression, as the 
20-pixel Sun was downsampled generally into 1-4 pixels and then 
compressed from 12 to 1 bit per pixel. For thumbnail images flagged 
with possible saturation, only non-thumbnail images could be used, 
and only after verification that saturation did not occur. 

The CCD response varies with wavelength and temperature 
(Bell et al., 2003). In general, blue filters show decreased response 
for increasing temperature, and near infrared filters show the 
opposite, as was seen in the pre-flight calibration for L8 and R8, 
respectively. Following Bell et al. (2003) but rearranging the terms 
slightly, we model the response as linear with temperature, such 
that F= F dns ■ R 0 ■ (1 + ■ (T+ 15 °C)), where F DNS is the measured 

total DN/s, R 0 is the responsivity at -15 °C, and Rj is the respon- 
sivity slope, T is the detector temperature in °C. A reference of 
-15 °C was chosen as the average temperature for images taken 
in the first Mars year of the mission. We define R 0 = 1 to maintain 
DN/s as the flux units of choice, but still account for the slope. 

The temperature slope for the CCD using the R8 filters was 
reasonably well determined in preflight calibration, but that using 
the L8 filters was not; its calculation is described here. The 
preflight calibration values relied on comparison with 430- and 
750-nm non-solar filters and an assumption regarding the red leak. 
Inspection of L8/R8 paired flux data shows that assumption to be 
wrong. If one assumes that x R = t l ■ (1 + e)— where e is a constant 
fractional difference in optical depth between the two filters— then 
the flux ratio for any given observation, FJF R , is approximately 
linear in temperature with a slope of Rjr-Rji and approximately 
linear in T sfan t (the slant-path optical depth, equivalent to x ■ rj) 
with a slope of e. We found that the average R ; difference was 
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0.00411 ±0.00002 and 0.00425 ± 0.00004 “(T 1 for Spirit and 
Opportunity respectively, while e averaged 0.0330 ± 0.0004 and 
0.0256 ± 0.0007. The implication of the larger temperature sensi- 
tivity of the CCD with L8 compared to Bell et al. (2003) is that 
the ratio of in- to out-of-band flux in L8 may be 6:4 rather than 
1:9, with the difference attributable to the relative blueness of 
the Sun compared to the calibration lamp and the fact that the 
in-band response was simply lost in the noise when imaging the 
calibration lamp. No trends or outliers were seen in the fit, suggest- 
ing that there were no detectable changes in e (for instance, due to 
the occasional presence of large ice particles). However, we calcu- 
late that e is only weakly size-dependent for dust with mean radius 
in the 1-5 pm range and a size distribution with variance 0.2, even 
without the red leak. 

3.3. Flat field 

Prior to summing pixels, a flat field image would normally be 
applied to account for high-frequency pixel-to-pixel variability 
and low-frequency variations across the detector caused by the 
optical system (intensity fall-off away from the optical axis, 
vignetting, out of focus dust on optical elements). However, accu- 
rate pre-flight solar-filter flat-field images could not be obtained. 
Certain data sets (Sol B/260, P2694; B/290, P2694; B/1276, 
P2694; B/1588, P2679; A/312, P2699; A/1295, P2691) imaged the 
Sun repeatedly, changing only the position of the Sun on the 
CCD, allowing a test of the flat-field strategy. Against several op- 
tions— nearest regular-filter flat-field, weighted-average of regu- 
lar-filter flat-fields, a 2-dimensional polynomial fit and cos 4 (<p) 
fall-off, where cp is the angle off the optical axis— defining the flat 
field as unity produced the lowest X 2 for sols <1200. With respect 
to high-frequency variations of order <1%, this is not problematic, 
as 300 pixels were summed when determining the solar flux. Expe- 
rience suggests that mid-frequency (one to several times the Sun’s 
apparent diameter) variations could be of order 1%, and may be 
caused by out of focus dust in the optics and other sources. For 
Opportunity, after the dust accumulation event described below, 
a second order polynomial fit in line and sample was used to cor- 
rect all observations for variations seen in P2694 and P2679. The 
corrections are up to 2% before the dust event and up to 23% after 
the dust event, increasing primarily to the right across the frame. 

3.4. Changing window dust 

The cycle of accumulation and removal of dust on the Pancam 
windows posed a problem in the determination of atmospheric 
opacity even for observations with the Sun in the same position. 
Interpretation of extinction from new window dust as new 
atmospheric dust was operationally significant, as the error would 
interfere with the assessment of the performance of the solar 
arrays. In addition, the retrieved window extinction came to 
exceed atmospheric extinction for parts of the mission. 

Accumulation of dust was recognized prior to the sol ~1 250 dust 
accumulation event. However, the accumulation event is illustra- 
tive. After the event— which occurred during the largest dust storm 
yet seen by the rovers— opacity was elevated. This was not diagnos- 
tic, as it could be considered routine when recovering from a plane- 
tary-scale dust cloud. However, diurnal variations— variations 
inversely correlated with airmass— were diagnostic. The process 
for identification and correction of dust events relies simply on a 
continuous reapplication of the Beer-Lambert Law. The transmis- 
sion of atmospheric dust decreases with airmass, while the trans- 
mission of dust on the window is independent of airmass, and 
mimics a lower value of F TOA . A 2-component model was used, where 
F = F T oa ■ exp( x window ^afmo'fl)i thus, ln(F) ln(Fyo/0 = —x wlrK i aw - 
- 4-atmo ■ >h Fig. 4 shows the resulting fit for every sol with at least 
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two Sun images taken over a range of >0.5 airmasses and at airmas- 
ses <5; error bars were calculated assuming the dominant error is 
natural variability in optical depth at the 5% level. A preliminary 
calibration was used to scale the data; for the final calibration, we as- 
serted r W jndow = 0 at landing and adjusted F TOA . That choice has no 
impact on the atmospheric optical depth retrieval. While no signif- 
icant degradation of the CCDs over the mission is expected, it would 
be part of the correction attributed to window dust. 

The history of dust on the window can be inferred from Fig. 4; 
so too can limitations on the method. High-frequency variations 
are typically within the noise. Although we associate the slower 
variations with dust build up, the higher-frequency variations 
could result from an error in the temperature calibration modu- 
lated by shifting patterns in time-of-sol sampling. We chose to 
conservatively interpret the data into a dust history, by modeling 
dust build up as a series of linear changes, with the slope changing 
at discrete times. For the model, the regimes were chosen initially 
by hand; then within each time range, a linear fit for window-dust 
as a function of time was determined. The fits were spliced to- 
gether, with their intersections defining the new regimes. Other 
functional forms, such as exponential decay, could be used, but 
we found none that provided a decrease in the residuals. 

All accumulation-regime boundaries were assumed to be con- 
tinuous in dust amount, but not in rate of accumulation. They 
are well-determined at the ~2% level, except for that near sol B/ 
1260. For that change, fits are unavailable for sols 1216-1280. Indi- 
vidual Sun images are available for most sols, with some gaps such 
as 1236-1242 at the height of the dust storm when the rover had 
insufficient power. Inspection of solar array data shows rapid dust 
accumulation over sols 1254-1268. If the Pancam dust accumula- 
tion occurred at that time, then the pre- and post-peak compari- 
sons of opacity and array data are in reasonable agreement. The 
retrieval was done using this assumption, but using an artificially 
large systematic error bar through this period The peak of the 
storm as seen in solar array data was sol 1237, when we estimate 
opacity at 4.8 from the array data, but when imaging was not 
commanded. 

3.5. Airmass 

Airmass is the ratio of path optical depth along an arbitrary line 
of sight to that in the zenith direction. Airmass is frequently 


approximated as sec(O), where 0 is the solar zenith angle; this rela- 
tion applies when the atmosphere thickness is small compared to 
the planetary radius. While we did not make that approximation, 
the distinction is immaterial for the vast majority of the data. We 
discuss here the minority of the data for which our assumptions 
regarding airmass matter. The extreme case is the moment of sun- 
set: infinite airmass for a plane-parallel atmosphere, but airmass of 
about 22 for spherical shell atmosphere with a scale-height of 
1 1 km and a radius of 3396 km. 

When the Sun is below an elevation angle of 30° (with respect 
to the gravitational normal vector), the airmass begins to be sensi- 
tive to the vertical distribution of the dust. Above that elevation, 
and the airmass for a plane-parallel atmosphere (e.g., all dust con- 
fined to a thin boundary layer such that the vertical scale is small 
compared to the planet’s radius) and that for an 1 1-km scale height 
differ by <1%. Below 10° elevation, the difference exceeds 10% and 
the airmass becomes somewhat sensitive to the specific choice of 
scale height. Other measures of vertical structure are possible: 
while these data are sensitive to vertical structure, they are not 
diagnostic. 

We determined an airmass relation to use for each rover by (1 ) 
finding all sols with valid data both with the Sun below 10° and 
above 30° elevation; (2) averaging all high-Sun (elevation >30°) 
data for each sol to get a characteristic mean opacity; (3) assuming 
the actual opacity when the Sun is low to be within ±5% (1-cr) of 
the mean; (4) taking the ratio of each observed slant opacity to 
the mean zenith opacity to be the airmass; (5) finding a single 
scale-height that best fits the observed distribution of airmasses. 
We chose a scale-height model in order to require only one new 
free parameter for vertical variation, due to the lack of uniqueness 
inherent in using airmass to infer vertical structure. The best fit 
scale-heights are 11.5 ±0.4 and 12.2 ±0.4 km for Spirit’s L8 and 
R8, respectively; and 19.3 ± 0.5 and 19.1 ± 0.5 km for Opportunity’s 
L8 and R8. For defining the elevation-airmass relation, we adopted 
model scale heights of 11.9 ± 0.4 and 19.2 ± 0.5 km, for Spirit and 
Opportunity respectively. We did not reduce the standard error 
values to account for averaging, as the sampling biases are similar 
for each filter. Implications are discussed in Section 4. The derived 
scale heights for Spirit are similar to typical gas scale heights near 
the surface. Those for Opportunity are surprisingly high, and likely 
arise from a real distribution that is different from a purely expo- 
nential model, as discussed in Section 4.5. 



1500 2000 

Opportunity sol number 


Fig. 4. Accumulation of window dust. Beer-Lambert law intercepts for individual sols are shown for each rover and filter, with all changes from sol 1 interpreted as 
accumulated opacity on the Pancam windows, which are physically separated for left and right eyes. Spirit L8 fits are at the top; Spirit R8, Opportunity L8 and R8 are 
successively offset by 1 unit optical depth. The abscissa is Opportunity sol number (Spirit sols are offset from Opportunity by -20.5). The ordinate is reversed to show 
accumulated optical depth as an absorption such that the Beer-Lambert intercept for a given sol is depressed: ln(F so i/F TOA ) = t. The window dust model is a 3-part line for 
Spirit (accumulation, rapid accumulation, removal) and an 8 part line for Opportunity. 
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3.6. Precision and accuracy 

Each individual measurement is a product of both a specific 
photon counting event and a relative calibration process. For 
nearly all data, counting statistics are irrelevant and systematic 
sources dominate the noise. The key relation can be written: 

Aslant = t -rj = ln(F TO A) - ln(F DN s) - ln(l + Ki ■ (T + 15°C)) 

— \n(Tfi at ) — i W i n d 0W . 

Therefore, when considering slant path opacity with a temper- 
ature correction <rl, 

& slant ~ [(ctroa/froA) 2 4- {lj dns / P dns ) 2 + (CRl • [T + 15°C]) 2 
+ (°flat/Tf,at) 2 + Window] V2 • 

And in general for the reported normal-incidence opacity, 

ffr = l[((y slant / t: slant? + (CT,,/^) 2 ] V2 . 

For most data, a z ~ ff S tant/<7. as the airmass uncertainty is negli- 
gible. For the other terms, counting noise was typically small: the 
ratio of noise to signal is approximately N _1/2 , where N is the total 
number of photons in the exposure, typically 6 x 10 5 . The read 
noise and uncertainty in the background subtraction were larger 
for low-signal images, and were explicitly considered in the error 
propagation. Uncertainty in the temperature term is a few times 
10~ 4 °C~ 1 [T + 15 °C], which, if all data were taken at -15 °C, would 
have no contribution. For data taken at CCD temperatures above or 
below the mean, this term contributes typical uncertainties of 
3-5 x 1CT 3 . For the flat field correction, Ofi at was assumed to be 
0.02, with Tfl at typically near unity, and this was initially the 
dominant instrumental source of uncertainty. 

Uncertainty due to non-photometric conditions (i.e., non-uni- 
form atmospheric dust) was the most difficult to quantify. We 
grouped the uncertainty associated with the top-of-atmosphere 
fit and the window dust model, as they have identical impact on 
a Beer-Lambert fit at any given time during the mission. Fig. 4 
shows error bars in individual fits, as determined from the uncer- 
tainties noted above and an assumption that natural dust-loading 
variability is 5%, which was determined to be consistent with the 
quality of the individual fits. Note that we omitted from these fits 
any sols during times of rapid sol-to-sol changes such as local dust 
storms or regional dust storm onset. While treating each measure- 
ment as independent would result in a lower estimate for the 
uncertainty in the fit, we took the largest systematic departures 
from the fit to imply the error associated with the top-of-atmo- 
sphere flux and window dust to be 0.03 (1-er) for each filter, except 
0.04 for Opportunity’s R8. In addition, we increased the uncertainty 
for times of rapid window-dust accumulation by the net change in 
10 sols; and we increase the uncertainty for sols A/1500-1800 and 
B/1216-1280 by 25% of the modeled window dust. These latter 
changes represent our estimate of a prudent systematic error to ac- 
count for times when the model is poorly constrained. 

The final source of uncertainty to note is that associated with 
airmass. Two factors are the uncertainty related to choice of model 
and the uncertainty relating to time. For the former, we assumed 
the scale-height model is an adequate description of the airmass 
variation with solar elevation, and we used the nominal scale- 
height uncertainties. For the latter, we adopted a time uncertainty 
of ±15 s based on the project requirement to maintain time knowl- 
edge (the translation from the spacecraft clock to spacecraft event 
time in UTC, from which the Sun’s position is calculated) of 30 s in 
the worst case. This translates to an elevation error of <0.1°. Air- 
mass errors are irrelevant when the Sun is above 1 5°, but dominate 
immediately before sunset when they approach 10%. 
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The resulting error represents systematic and random effects as 
discussed above. In most cases, an absolute error of ~0.03 /i) 
(~0.04/i; for Opportunity’s R8) can be assumed to be system- 
atic-resisting attempts to use averaging to reduce error— while 
the (quadratic) remainder can be assumed to be random. Errors 
in high-airmass (>10) optical depths are dominated by systematic 
effects. Averaging to reduce the error is not recommended due to 
the identified systematic terms. A conservative approach would 
be to make comparisons across and within sols with consideration 
for their (quadratic) cumulative error. 

3.7. Planetary Data System archive 

The data are archived with the Planetary Data System (Lemmon, 
M., MER Mars Pancam Atmospheric Opacity RDR VI. 0, NASA Plan- 
etary Data System, MER-M-PANCAM-5-ATMOS-OPACITY-Vl.O, 
2004). There are two files for each combination of rover and cam- 
era, a label (*.LBL) and tabulated data (*.TAB). An example filename 
is “2TAU_440_2209_201 1021 8A.TAB", where “2” signifies MER-A, 
Spirit (“1” is MER-B, Opportunity); “440” is the nominal wave- 
length of the L8 filter (“880" is the value for R8); “2209” is the last 
sol appearing in the table; “20110218” is a string signifying the ta- 
ble creation date (YYYYMMDD). As is customary, the labels de- 
scribe the contents of the tables. The tables themselves include 
columns for EDR product identifier, L s , Mars-Sun distance, decimal 
sols since midnight before landing (sol number-1 +LTST/24— one 
rounds up to get project sol number), airmass, flux in corrected 
DN/ms, optical depth, and 1-tr error in optical depth. 

The results discussed here refer to the calibration associated 
with the time stamps “20110218” for Spirit and “20131028” for 
Opportunity. Future deliveries will be updated according to contin- 
ued application of the techniques described here. 

4. Discussion 

4.1. Seasonal variations of optical depth 

Fig. 1 shows the optical depth record of the mission. Values in 
both filters are nearly the same: values discussed in this section 
are from R8 (880-nm) unless specified otherwise. Due to the 
well-understood spectral sampling, the R8 opacity is recom- 
mended as the value of reference. The rovers’ initial 90-sol mission 
included the southern autumnal equinox (L s = 0°) that began MY 
27, during a time of moderate but declining optical depths 
(Lemmon et al„ 2004) following a regional dust storm that was 
dissipating by the landing time (Smith, 2008). The TES climatology 
(Smith, 2004) would suggest that optical depths would decline into 
aphelion season for both rovers, rise slightly near L s = 135°, and 
become dusty during perihelion season, with dust increases 
around L s = 180°, L s = 220°, L s = 320°. Spirit optical depths declined 
from near 0.90 to below ~0.3 by sol 155 (L s = 45°), and remained 
similarly low until about sol 350 (L s =135°). During this time, 
Opportunity optical depths declined from 0.95 to below ~0.5. 
The sudden increase of optical depth over L s = 135-137° was an 
operational surprise. It occurred in the last week of 2004, during 
a time of reduced staffing and was seen as the first major reduction 
in solar power of the mission: an optical depth 1.24 dust storm 
peaking over sols B/329-332. Both rovers experienced local dust 
storms with optical depths of 1-1.5 and generally elevated optical 
depths from L s = 135° to 200° (to about sol 500 for Spirit). 

For context, Mars Global Surveyor (MGS) Mars Orbiter Camera 
(MOC) views of some of these storms are shown in Fig. 5. The 
MGS spacecraft reached its designated 380-km altitude, 1400 LMST 
equator crossing (south-to-north), circular polar orbit in March 
1999 and operated through November 2006. For that time, MOC 
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provided nearly continuous daily global mapping of the martian 
surface at 7.5 km pixel -1 in two band passes, red (600-650 nm) 
and blue (400-450 nm) (Malin et al., 2010). In the spirit of mutual 
collaboration between missions, MOC and later the Mars Recon- 
naissance Orbiter (MRO) Mars Color Imager (MARC1) provided al- 
most daily regional weather updates and forecasts for the two 
MER sites during dusty seasons. 

From late December 2004 (L s =131°) through March 2005 
(L s = 185°) there were several significant dust events, which either 
passed through or near the Meridiani and Gusev Crater rover sites. 
The earliest was the most extensive of these dust events, observed 
as the MGS spacecraft came out of contingency mode on 27 
December 2004, L s = 135° (Cantor, 2007). Prior to that time, weath- 
er conditions at southern low-latitudes were normal for mid-win- 
ter. The aphelion water-ice cloud belt was at its maximum extent, 
orographic water-ice clouds were observed over the Tharsis volca- 
noes, and an occasional local dust storm in noted Syria. By 
L s = 135.5°, dust storm activity extended from 0° to 40°S, 300° to 
90°W with most of the activity occurring in Noachis and Valles 
Marineris, as well as locally in Margaritifer. Storm activity ex- 
tended into Meridiani by L s = 137.0°; see Fig. 5a. The origin of the 
dust event was not known exactly, due to the absence of MOC 
observations during the onset of storm activity. But it was clear 
from the first observations, that the dust event was not a single 
storm but consisted of at least one regional storm in Noachis with 
local storms in Valles Marineris, Margaritifier, and Meridiani. The 
storm activity was extensive enough to generate a diffuse dust 
haze that encircled the southern hemisphere within a two week 
period, raising the visible background opacity above nominal levels 
for this time of year, to 0.7 at both sites. Cantor (2007) noted that 
this event was seasonally the earliest planet-encircling dust veil 
observed to date. Dust-lifting activity abated by L s = 138.5°, only 
several sols after it was first observed. However, the background 
opacity in the southern hemisphere remained elevated for over a 
month, which resulted in water ice clouds associated with the aph- 
elion cloud belt and the shield volcanoes to disappear for that per- 
iod of time. 

Within the 70 sols after the planet-encircling dust veil, a series 
of dust storms passed over the Spirit site. The first was observed by 
MGS on the 19th of January 2005 ( L s = 146.5°); see Fig 5b. Initially 
the local storm formed in Promethei along the south polar cap 
edge. The storm moved east-northeast at ~13 m s -1 over a 2-sol 
period before crossing over Gusev crater. The storm quickly abated 
the following sol, but a diffuse dust haze remained over the site for 
several additional sols. The next storm fronts to reach Gusev were 


observed from 6 to 10 March (L s = 170.6-173.0°, sols 417-421). 
That activity consisted of two storms, a single regional storm that 
passed through the Spirit site on 7 March and a large local storm 
that came through two sols later. Both storms began along the edge 
of the seasonal receding south polar cap edge. The regional storm 
began as a large local storm in Promethei on L s = 168° and over a 
2-sol period moved east at an average speed of about 22 ms -1 , 
merging with local storms along the cap edge and dramatically 
increasing in area. The storm reached its maximum areal extent 
(~6.4 x 10 6 km 2 ) at L s =169°. At that same time, the leading 
edge of the storm was in Sirenum, far to the southeast, while the 
western trailing edge extended south-southwest of Gusev. The 
western, trailing edge of the storm moved north over the Spirit site 
at L s = 171°. The storm was gone as quickly as it came, moving off 
to the east by the following sol. The second storm, developed just 
to the south of Gusev in Cimmeria. The storm moved northward at 
about ~13.5 m s -1 , completely enveloping the Spirit site at 
L s = 173°; see Fig. 5c. The storm, like the previous local storms, 
had dissipated by the following sol, but it left behind a diffuse dust 
haze that persisted for several additional sols. 

Near L s = 215-220°, a new period of rising opacity was seen first 
at the Spirit site, but briefly peaked at optical depth 1.98 at the 
Opportunity site on B/492, as cross-equatorial storms on the Acida- 
lia storm track created a second planet-encircling haze (Cantor, 
2007). Optical depths fell gradually until L s = 307° before the next 
wave of cross-equatorial storms increased opacity to 1.76 on sol B/ 
630 (L s = 316.5°) and eventually encircled the planet. A Hubble 
Space Telescope (HST) Advanced Camera for Surveys (ACS (Sirianni 
et al., 2005)) image (Fig. 6) shows a 1500-km long regional dust 
storm in Meridiani on sol 626 (L s = 314.2°, 28 October 2005) just 
before the storm peaked at the Opportunity site. After this time, 
optical depths at both sites gradually declined while the mission 
entered its second Mars year. By L s = 0°, the observations essen- 
tially matched those of the initial 90-sol mission. 

Subsequent Mars years displayed broadly similar seasonal 
behavior, with the major exception being the globally extensive 
dust storms around L s = 270° of MY 28. Fig. 7 shows optical depths 
for all years in their seasonal context. The optical depths generally 
match orbital climatology from imaging and infrared opacity, 
which show low optical depths through southern winter, receding 
cap-edge dust storms leading to an increase in dust after L s = 135°, 
and regional storms around L s = 220° and 315° (Smith, 2004, 2009; 
Wang, 2007). 

The period from L s = 0-135° was substantially the same in 
each year. Optical depths at the Spirit site show a fall-off over 



Fig. 5. Local dust storm activity at the two MER sites. (A) In Meridiani Planum on sol B/330 (28 December 2004, L s = 135°) the optical depth was near its peak of 1.2. (B) In 
Gusev Crater on sol A/372 (19 January 2005, L s = 147°) the optical depth was 1.1, and peaked at 1.4 the next sol, with slightly elevated optical depths after the 2-sol event. (C) 
In Gusev crater on sol A/420 (9 March 2005, L s = 173°) the optical depth peaked at 1.4 in the afternoon, and varied from 0.7 to 1.5 within a 10-sol interval. Maps are composed 
of MOC daily global mapping images that have been cylindrically map projected at 6 km pixel -1 and mosaicked together. The color has been stretched to enhance the contrast 
between the storms and the martian surface. The Meridiani frame spans 27.2°S-22.8°N, 154.0°-214.0°W and the two Gusev frames span 40.0°S-10.0°N, 155.0°-215.0°W. The 
small circle near the center of each image frame corresponds to the location of the rovers. (For interpretation of the references to color in this figure legend, the reader is 
referred to the web version of this article.) 
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Fig. 6. Regional dust storm activity entering the Meridiani region. (A) A global view was acquired by the Hubble Space Telescope’s Advanced Camera for Surveys, using 250-, 
502-, and 658-nrn filters, on 28 October 2005 (L s = 314°, sol B/626). The area of the 1500-km long storm is shown (B) in a 2500-km wide excerpt of the global view and (C) in a 
comparable Wide Field Planetary Camera 2 image acquired 26 June 2001 with storm-free conditions. The circle shows the Opportunity site, where the optical depth was 1.4 
and peaked at 1.8, 4 sols later. [Image credit: NASA, ESA, The Hubble Heritage Team (STScI/AURA), J. Bell (Cornell University) and M. Wolff (Space Science Institute).] 


L s = 10-50°, followed by a slight rise through L s = 70, a decline to 
L s = 85°, a rise to I s = 105°, and a decline to L s = 135°. The features 
are low-amplitude, but repeatable, and suggest an influx of dust 
around 15-20° before and after northern summer solstice. The 
most notable departures occur during MY 29— few calibration 
observations were able to be acquired during this time, and the 
systematic uncertainties are at their largest. It is reasonable, but 
not necessary, to suppose that the calibration is off and MY 29 
should be offset slightly lower to mimic the other years. At the 
Opportunity site, the first year’s pattern repeats, including an 
inflection in the slope of decline at L s = 75°. During this time, 
several individual data points he above the trend as seen with both 
filters. These are likely due to discrete clouds, as discussed in 
Section 4.3. 

The period from L s = 135-215° was similar to the first year. The 
optical depths did not repeat as they did in the winter, but the per- 
iod was dominated by generally increased opacity. The period did 


not start with as sudden a rise as in the first year, but the rovers 
typically experienced ~1 small storm of a few sols in duration 
and peak optical depths of over 1 to about 2. MY 28 is notable in 
showing a more gradual increase in Meridiani than the other Mars 
years. The dust load away from local storms tended to vary by tens 
of percent across years, and the timing and number of local storms 
(or brief opacity spikes) varied. 

For L s = 215-315°, MY 29 is almost identical to MY 27, with re- 
gional dust storms increasing opacity at both sites early in the per- 
iod and then a general decline. MY 30 and 31 are similar at the 
Opportunity site, with a later onset of diffuse opacity from regional 
storms and smaller opacity increases. For years other than MY 28, 
the optical depths from L s = 240-305° are remarkably consistent, 
especially for Opportunity (Fig. 8). 

MY 28 is distinguished by a somewhat steady and relatively low 
optical depth around L s = 215-265°, followed by a rapid increase 
beginning at L s = 265° for Opportunity and 269° for Spirit (26 June 
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Fig. 7. 880-nm optical depths as a function of season. Shown are results from (A) 
Spirit and (B) Opportunity for Mars years 26 (cyan, L s > 330 only), 27 (black), 28 
(red), 29 (green), 30 (blue), and 31 (magenta). The 2007 dust storm is cut off to 
show detail at other times. (For interpretation of the references to color in this 
figure legend, the reader is referred to the web version of this article.) 






Ls 

Fig. 8. 880-nm optical depths for southern summer and autumn. Conventions are 
the same as Fig. 7, except the logarithmic axis shows the full range of optical depths. 
Solid maroon lines illustrate (in A) a 2.3%/sol decay over L s = 290-330°, and (in B) a 
2.4%/sol decay over L s = 285-315° and a 1.3%/sol decay over L s = 315-335°. 


and 2 July 2007). This storm led to optical depths up to 4.3 at the 
Spirit site, where the optical depth was about 3 or more for a 30-sol 
period and rover operations were dramatically slowed by lack of 
available solar power. At the Opportunity site, the observed peak 
optical depth of 4.6 occurred on sol 1235 (L s = 277, 16 July 2007), 
after which the rover was commanded to hibernate to slow the 
drain of the batteries. The next observations, on sols 1243 and 
1246, showed optical depths down to 3.9 and then back up, to 
eventually reach 4.5 and then slowly decline. Analysis of solar ar- 
ray data suggests the peak sol-averaged optical depth was proba- 
bly near 5, comparable to estimates of the first VL1 storm, but 
short of the second, t = 9 storm (Pollack et al., 1979). MRO MARCI 
(Bell et al., 2009) images were used to track this event globally. 
Fig. 9 shows the Opportunity site as a series of individual storms 
pass near the site and the general dustiness of the atmosphere in- 
creases. After L s = 295°, both sites show a similar, steady decrease 
of opacity (see Fig. 8). The initial rate of decrease is 2.3-2.4%/sol, 
before the decrease slows and the optical depths become similar 
to those from other years. Smaller dust storms typically had higher 
peaks rates of decline maintained for shorter durations, indicating 
advection as a likely means of removal. Global dust storms re- 
ported by Viking initiated near L s = 205° and reached optical depth 
2. 7-3.2 at the VL1 site; and L s = 270°, reaching optical depth 3.7-9 


attheVLl site (Pollack et al., 1979). Our decay timescales of 43 sols 
for each rover and then 77 sols for Opportunity are similar to the 
decay rates following the VL1 second global storm of 51 and then 
119 sols (Pollack et al., 1979). The lower rate, which was observed 
for both VL1 storms but not for Spirit, is likely due to renewed dust 
lifting in the Acidalia cross-equatorial storm track. 

As observed from orbiters, the period from L s = 315-360° typi- 
cally produces a final wave of cross-equatorial dust storms (Cantor 
et al., 2001; Wang et al„ 2003; Wang, 2007; Cantor, 2007), result- 
ing in elevated dust optical depths followed by settling at the rover 
sites. This activity, as seen at the rover sites, was sporadic, with de- 
cay of the global event at both sites in MY 28, regional storms at 
the Spirit site in MY 29 and the Opportunity site in MY 27, 29, 
and 31 (and prior to landing), and general lifting otherwise. 

The overall picture is qualitatively the same as previous reports. 
Smith (2004) presents annual and zonal averaged TES-derived 
optical depths at 9 pm for MY 24-26. The seasonal trends are the 
same as those seen by the rovers. A gradual decline in dust optical 
depth to L s = 140° is followed by elevated backgrounds, then fur- 
ther enhanced storms (and the 2007 global dust event) after 
around L s = 220°, and then a third wave of lifting near L s = 340°. 
Smith (2009) presents a similar data set for MY 26-29 from THE- 
MIS, substantially overlapping the rover mission. These show the 
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Fig. 9. The 2007 global dust event at the Opportunity site. Shown are excerpts, spanning from 60°S-30°N, 330°-60°W, from Mars Reconnaissance Orbiter Mars Color Imager 
daily mosaics over 21 June to 18 July 2007 (sols B/l 2 1 1—1237, L s = 262-279°), at a (summed-8) resolution of 7.5 km/pixel. The white circle in each image shows the 
Opportunity site; black areas are gaps in data. During 21 -26 June (A-F), atmospheric optical depth increased from 0.9 to 1.4 and became variable on time scales <1 sol. (G)The 
atmospheric optical depth was 2.8 on sol 1223, 3 July; (FI) 3. 5-3. 9 (variable) on sol 1233, 14 July; and (I) unmeasured on sol 1237, 18 July. During sols 1236-1241, the rover 
was shut down, but solar panel output suggests a peak optical depth near 5 on sol 1237. (For interpretation of the references to color in this figure legend, the reader is 
referred to the web version of this article.) 


same pattern. In both cases, the 9-pm data are about half the Spirit 
observations at the same latitude. The differences between the 
data sets come mostly from the visible vs. infrared properties of 
the dust and the presence of water ice at the Meridiani site in 
the winter (both discussed below) and slightly from meridional 
variations. 

Away from the globally extensive dust event, the opacity de- 
clines around L s = 0, 180, and 270° occur when lifting and transport 
may supply or remove some dust. We considered every 10-sol per- 
iod in the mission in 5-sol increments, measured the goodness of 
fit of an exponential change model and the associated rate of 
change. To investigate periods of steady decline, we exclude peri- 
ods for which the data were poorly fit (such as the passage of local 


storms) and periods that had dust increasing or stable (slope 
consistent with 0 or less at the 90% confidence level). We found 
remarkable stability for each rover outside of the storms men- 
tioned above. For Opportunity, the average change during periods 
of steady decline was 0.52 ± 0.05% per sol; for Spirit the decline 
was 0.77 ± 0.06% per sol. The rates and their implied removal 
timescales of 130-200 sols are consistent with initial reports of 
0.6-0.7% per sol by Lemmon et al. (2004) for the spring period. The 
rate is also consistent with the settling rate determined for 
Pathfinder of 0.28%/sol by surface coverage (Landis and Jenkins, 2000) 
during a period of optical depth ~0.5 (Smith and Lemmon, 1999), 
implying an atmospheric loss to sedimentation of 0.56%/sol. Thus 
sedimentation is likely a significant component of dust loss 


Please cite this article in press as: Lemmon, M.T., et al. Dust aerosol, clouds, and the atmospheric optical depth record over 5 Mars years of the Mars 
Exploration Rover mission. Icarus (2014), http://dx.doi.Org/10.1016/j.icarus.2014.03.029 






ARTICLE IN PRESS 


M.T. Lemmon et al./Icarus xxx (2014) xxx-xxx 


12 

through much of the year. Given the timescales above, areas that 
are net sinks for atmospheric dust could gain up to 3-5 times the 
mean atmospheric dust load each Mars year. 

4.2. Seasonal variations in aerosols 

In principle, varying color ratios can be indicative of aerosol size 
changes, and color ratio (or Angstrom exponent) has some diag- 
nostic value. We find that (1) the ratio of the L8 to R8 opacity does 
not show significant variation during the first part of the mission; 
and (2) variability late in the mission is unmeasurable due to 
uncertainties in the correction of variable dust on the windows. 
Section 3.2 reports the mean ratio of (R8 - L8)/R8 opacity as 
0.0330 ± 0.0004 and 0.0256 ± 0.0007. With the nominal factor of 
two wavelength difference, the ratios correspond to an Angstrom 
exponent near zero, typical of large particles. Dust aerosol size dis- 
tributions have been reported with sizes around 1.5-1.65 pm (ra- 
dius of equivalent volume sphere) and variances around 0.2-0.5 
for Viking, Pathfinder and MER (Pollack et al., 1995; Tomasko 
et al., 1999; Lemmon et al., 2004). Such aerosols have an extinction 
maximum near or longward of 1000 nm, and have about 7-11% 
less opacity at 440 nm than at 880 nm. Given the red leak for the 
L8 filter, we cannot exclude such aerosols. 

For the first part of the mission we have measurements of 9-pm 
dust opacity from Mini Thermal Emission Spectrometer (Mini-TES) 
data, as reported by Smith et al. (2004, 2006). The record was con- 
tinued through the middle of MY 28 before dust on the Mini-TES 
optics rendered the data useless. Fig. 10 shows the Mini-TES data 
with the 880-nm data for both rovers— note the infrared opacity 
scale is half the visible scale to generally overlay the data. All data 
are shown except three measurements near optical depth two on 
sol A/1243. These were two times the 880-nm measurements at 
the time, but are subject to increased uncertainty due to the dusty 
optics. The infrared/880-nm ratio is shown in Fig. 10(C). Conven- 
tion is to report the inverse as a “visible/IR” ratio. However, due 
to the distribution of error bars, the figure is more interpretable 
showing “IR/visible”. There are clear trends: after landing, the ratio 
falls for both rovers; Opportunity shows more complex behavior 
during the first northern summer; the ratio increases significantly 
with the onset of dust storms; the ratio tends to decline after dust 
storms and generally increases again for new storms; the ratio 
then falls for the next northern summer and rises with the onset 
of the second dusty season. Interestingly, the L s = 220° dust storm 
did not change the ratio at the Spirit site. Note that the ratio of sol- 
average optical depths is shown, so the few sols during which there 
was notable variability may have incorrect ratios. 

The comparison of these optical depths can be informative as to 
dust aerosol size (Clancy et al., 2003; Wolff et al., 2006). For certain 
assumptions about IR optical properties of dust, Wolff et al. (2006) 
derived the (inverse) ratio as a function of mean size and variance 
of the size distribution. Based on their Fig. 13 and assuming a var- 
iance of 0.5, in Fig. 10(C) we show the ratio for mean radii of 0.7, 
1.4, and 2.1 pm for both rovers. For different assumptions, those 
sizes would change, but the approximate relation would remain. 
The Opportunity data around L s = 30-130° are somewhat problem- 
atic: the 880-nm value is a total optical depth; the Mini-TES value 
is specific to dust. Mini-TES did not detect water ice, but orbital 
data show ice at that time, and the Mini-TES data do not rule ice 
out at high altitude (Smith et al., 2006). So, the ratio may be de- 
pressed due to ice hazes. It is possible that the increase in the 
Mini-TES /880-nm ratio during this period may indicate a sensitiv- 
ity to ice-coated dust. Generally, dust sizes fall as low as ~1 pm in 
northern summer and increase to over 2 pm with the onset of a 
dust storm (given the caveat above). This pattern is similar to pat- 
terns reported by Clancy et al. (2003) using orbital infrared spec- 
troscopy, and generally consistent with dust sizes of 1.2-1. 4 pm 


during the northern summer at the northerly Phoenix site 

(Komguem et al., 2013). 

4.3. Clouds 

The Spirit and Opportunity optical depth records are quite sim- 
ilar. One systematic exception is the elevated opacity during north- 
ern summer at the more northerly latitude. Additionally, during 
the L s = 0-135° time, Opportunity optical depths tend to steadily 
decline, with infrequent outliers (in both filters) of magnitude 
~0.1 above the trend. While Mini-TES did not directly detect ice 
aerosols, possibly due to their location at cold, high altitudes, the 
TES climatology and contemporaneous orbital observations do 
show near-equatorial ice aerosols (Smith, 2004, 2009). The TES cli- 
matology shows slightly more dust at the Opportunity site com- 
pared to the Spirit site, but also shows much more ice, during 
the L s = 0-150° timeframe, with a sudden decrease at L s = 135° 
during the first year. For each rover, this is the only time of year 
when appreciable ice-aerosol amounts are expected. While separa- 
tion of ice and dust was a motivation for inclusion of Pancam’s two 
solar filters (Bell et al., 2003), the sensitivity is only to particles of 
very different sizes, and is significantly reduced due to the L8 red 
leak. 

During the first two years (after sol 150) regular sky monitoring 
was done with the rovers’ Navcams, taking advantage of the 
45 x 45° field of view (Maki et al., 2003). Throughout the missions, 
regular drive-direction and other panoramic imagery frequently 
showed the sky. Under typical aphelion-season conditions, we ex- 
pect line-of-sight opacity variations in excess of a few percent to be 
detectable, although the limit does increase for dustier skies. Prior 
to sol B/100, no discrete clouds were noted. During sols B/101-109 
(L s ~ 30°) and again in the 140 s (L s ~ 50°), clouds were noted in 
several images, leading to the regular campaigns (Fig. 11). During 
the first Mars year, cirriform clouds were seen occasionally until 
sol B/310 (over L s = 50-126°), especially around L s ~ 55° and 
115°. During the second year, cirriform clouds were seen over 
L s = 23-50° and L s = 110-131°. During the third year, sampling 
was more sparse, but such clouds were seen near L s = 20° and over 
L s = 109-136°, with a real gap in between. Additional clouds were 
seen near L s = 30° and L s = 63° of the next two years (successively), 
but sampling had been curtailed due to operational constraints. 
Such clouds are likely the reason for the small spikes in optical 
depth. Obviously-cirriform clouds were not observed at the Spirit 
site, although some of the infrequent wave structures were likely 
due to condensates. 

The images are not diagnostic as to the cloud composition. 
Water ice is likely, given the aphelion cloud belt over the Opportu- 
nity site. However, mesospheric carbon dioxide ice clouds have 
been seen near the Meridiani area, peaking near L s = 30° and 
150° (Clancy et al., 2007). It is possible that some of the early- 
season clouds were C0 2 ice. Fig. 1 1 shows some of the variety of 
clouds, including wave and more complicated structures, not un- 
like terrestrial cirrus. 

The cloud detections were episodic, with many over a few sols 
and then many sols in between them. Yet the elevated opacity was 
consistent. Thus, the general opacity at the Opportunity site during 
northern summer was either dust, or dust with a somewhat uni- 
form ice haze that was not detectable as discrete clouds. No dis- 
crete ice clouds were seen at either site outside of L s = 20-136°. 

4.4. Seasonal variations in insolation and dust devils 

Dust devils are seen across Mars in orbital imagery (Malin and 
Edgett, 2001 ; Cantor et al., 2006), and they have been reported and 
used to estimate local dust lifting rates at the Mars Pathfinder, 
Spirit, and Phoenix landing sites (Metzger et al., 1999; Ferri et al., 
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Fig. 10. Visible to infrared opacity comparison. (A) The 0.88-pm (R8) opacity (black line) is compared with dust opacity observed by MiniTES at 9-pm (red diamonds, see 
Smith et al. (2006)). Uncertainties, not shown, are typically small compared to the symbols. (B) The same is shown for Opportunity, in blue. (C) The ratio of MiniTES to R8 
opacity is shown; every 20th error bar is plotted. Opportunity values are offset by 0.5. The dotted lines correspond to aerosol size distributions with variance 0.5 and mean 
radius of 0.7, 1 .4, and 2.1 pm, from bottom (see text). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


2003; Greeley et al., 2006a, b, 2010; Ellehoj et al., 2010). We note 
that in Opportunity images, dust devils do appear and dust is seen 
being lifted by wind gusts in craters; however, such occurrences 
are rare compared to the Spirit site, presumably due to the lower 
availability of surface dust and the generally smoother surface. 

Heating of the surface and atmosphere by sunlight is governed 
by both orbital parameters and atmospheric extinction. Mars has 
an eccentric orbit, with aphelion at 1.67 AU and perihelion at 
1.38 AU. Thus, globally averaged insolation at the top-of-atmo- 
sphere (TOA) varies by nearly 50% between these two orbital 
positions. The Opportunity site is nearly equatorial, so seasonal 
effects are dominated by the Sun being lower in the sky near 
the solstices. The Spirit site is still “tropical” in the sense of having 
the sub-solar latitude sweep over its latitude before and after 
summer solstice. However, it is sufficiently far south that the 
Sun is much lower in the sky in winter. The two effects are 
roughly in phase, with perihelion at L s = 251°, roughly 35 sols 
before summer solstice. 

We have modeled the solar fluxes at TOA, at the surface and ab- 
sorbed into the atmosphere (Fig. 12). The TOA calculations are 
straightforward geometric calculations, taking into account the 
varying distance to the Sun, length of daytime, and elevation of 
the Sun throughout the day at each site. Atmospheric calculations 
were done with the Tomasko et al. (1999) dust and radiative trans- 
fer model, supplemented by expanded dust parameters provided 
Johnson et al. (2003), every 15 min for the duration of the mission. 
The sol-averaged optical depth value was used in each calculation. 
Wavelengths outside of 440-1000 nm were extrapolated from the 
above parameters and contributed a minority of the flux. Atmo- 
spheric gas absorptions were ignored, and we consider only solar 


radiation, rather than thermal radiation. These models are similar 
to project-internal models used to monitor the solar panels and 
to estimate power loss due to dust deposition. One distinction is 
that this paper does not consider the effect of rover tilt and instead 
focuses on insolation over a level surface. 

Surface insolation is shown as flux incident on the surface, 
roughly 80% of which is absorbed at the Spirit site, 86% at the 
Opportunity site (Bell et al., 2008). Insolation follows the seasonal 
TOA trends, but is modulated by dust in the atmosphere. Even dur- 
ing the optical depth >4 dust storm, surface insolation only falls by 
a factor of ~2-3 (rather than ~100 for direct sunlight) due to scat- 
tered light reaching the surface. The figure shows total atmo- 
spheric column absorption, which is roughly correlated with 
optical depth. When light interacts with dust, the most likely out- 
come is forward scattering, so the proportionality of absorption to 
abundance holds to optical depths ~2; at higher optical depths, 
absorption increases less rapidly. For low optical depths, similar 
patterns would be seen at all levels of the atmosphere. For the 
highest optical depths, absorption is shifted higher in the atmo- 
sphere; however, this is not an especially strong effect due to the 
forward scattering. High levels of dust could be expected to play 
a stronger role in cooling the bottom of the atmosphere in the 
infrared, through increased radiative coupling to a cooled surface 
and to space. Note that L s = 135° dust storms can reduce the solar 
flux to mid-winter levels at each site. The global dust event re- 
duced flux at the Opportunity site to levels far below those in win- 
ter; Spirit did not see such a low value due to the near-solstice 
timing and a more southerly latitude. 

Greeley et al. (2010) discuss 3 martian “dust devil seasons” 
that occurred during the 3 southern summers of the Spirit 
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Fig. 11. Opportunity Navcam images of clouds acquired on sol (A) 269, L s = 106°; (B) 290, 116°; (C) 291, 117°; (D) 758, 24°; (E) 763, 27°; (F) 949, 112°; (G) 950, 112°; (H) 1647, 
126°; and (1) 2847, 63°. All images have a linear contrast-stretch; clouds are typically near 15% contrast at their peak, except for (B) 30% and (1) 5%. 


mission. Fig. 12(A) shows their estimates of dust devil density. We 
have added error bars and upper limits based only on the count- 
ing statistics reported in Greeley et al. (2010, Fig. 3), with no 
adjustment for time of day effects or any bias resulting from the 
detection method. We note that at the site, and for dust devils 
visible to the rover’s Navcam during the warm time of day 
sampled in the Greeley et al. data, dust devil number density 
seems to increase exponentially with surface insolation (a similar 
fit could have been achieved for noon insolation or for the 
surface-atmosphere difference) such that the dust devil density 
is lCT 3 km~ 2 at 116Wm~ 2 with an e-folding scale of 18 W itT 2 . 
Other forms, such as a linear fit with a threshold near 
140 W m~ 2 , are not ruled out due to the scarcity of reported data 
away from “dust devil season”. The exponential relationship 
between dust devil number density and mean surface insolation 
fits the seasonal change, the effect of the MY 28 dust storms, 
and the fall off of dust devils with the pre-solstice ( L s ~ 220°) dust 
storms of MY 27 and 29 and the final major storm of the mission 
(L s ~ 320°). Unlike a threshold model, the exponential is consis- 
tent with the formation of a new dust devil track near the rover 
at L s =17±5° during the first Mars year, far from “dust devil 
season” [M.C. Malin et al„ Wheel tracks from landing site to hills, 
NASA’s Planetary Photojournal (http://photojournal.jpl.nasa.gov/), 
PIA07192, 3 January 2005], The specific relationship cannot be 
extrapolated across Mars (or even to the Opportunity site), as dust 
supply and local boundary layer meteorology are also integral 
factors in dust devil formation. 


4.5. Vertical structure 

To process data taken with the Sun at a low elevation, we in- 
voked an atmospheric model with exponential fall-off of dust with 
altitude, such that the vertical distribution is fully expressed as a 
scale height in Section 3.5. This convenient 1-parameter model is 
adequate with respect to the data, but it provides limited con- 
straints on the vertical behavior of the dust. More complex dust 
structures (Whiteway et al., 2009; Heavens et al., 2011) are not ru- 
led out. The Spirit data are characterized by a scale height compa- 
rable to what one would estimate for a martian C0 2 atmosphere 
with a mean temperature of 210 K. The derived Opportunity scale 
height is significantly larger than the gas scale height. Such a scale 
height is unlikely, as it would imply a dust source high in the atmo- 
sphere and sink at the surface. However, any model that fit the 
data would share the characteristic that dust is preferentially 
farther away from the surface than would be expected for a well- 
mixed atmosphere. For example, a possible alternative is that, 
while the Gusev atmosphere has dust well-mixed for the bottom 
2-3 scale heights (the region the model is sensitive to), the Merid- 
iani atmosphere is top-heavy, perhaps reflecting the high-altitude 
tropical dust maximum at 20-30 km (Heavens et al., 2011). 
Alternatively, a relative depletion of dust in the boundary layer is 
consistent with the relatively lower dustiness of the Opportunity 
site relative to the Spirit site. In any case, while one can conclude 
the dust is relatively higher at the Opportunity site, the form of 
model applied is not uniquely constrained. 
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3. Ice clouds contributed to the observed opacity at the Opportu- 
nity site. Water ice hazes may have been present. Clouds were 
seen over L s = 20-136°, with peak activity near L s = 50° and 
115°. Ice clouds and probably hazes were not a significant part 
of the opacity at the Spirit site. 

4. The shortwave optical depth of the atmosphere drives impor- 
tant changes to surface insolation and to atmospheric heat 
input. Dust devils are common at the Spirit site when insolation 
is high, but are rarer when insolation is seasonally low or low 
due to dust storms. While the occurrence of dust devils in an 
area responds to many factors, the relative abundance across 
seasons correlates with surface insolation. 

5. Data from the Spirit site are consistent with the dust being well 
mixed with the atmosphere in the bottom 10-20 km. The 
Opportunity data require relatively more dust at higher 
altitudes, perhaps reflecting the high-altitude tropical dust 
maximum at 20-30 km. 



Fig. 12. Insolation at the rover sites for (A) Spirit and (B) Opportunity. The 
continuous curves show modeled top-of-atmosphere, sol-average insolation 
(upper, black), direct plus diffuse surface insolation (middle, blue) and atmospheric 
absorption of sunlight (lower, red). For Spirit, symbols (green) show dust devil 
number density (right axis) reported by Greeley et al. (2010) with T symbols 
indicating upper limits. (For interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this article.) 

5. Conclusions 

1. A precise optical depth record for the Mars Exploration Rover 
sites has been created, using direct solar extinction measure- 
ments, and archived. The data set spans more than 3 Mars years 
for Spirit and more than 5 Mars years for Opportunity, including 
periods of local, regional and globally extensive dust storms. 
The data set shows that, for each site, the L s = 0-135° period 
is consistent year-to-year, with relatively lower dust amount;, 
the L s = 135-215° period has variable local to regional storms; 
the L s = 215-360° includes a global dust storm, but otherwise 
shows a similar set of dust lifted by regional storms. The 
880-nm (Pancam R8 filter) optical depth should be used for 
comparisons among data sets and model inputs, as the L8 
(440-nm) filter has a significant 750-nm contribution. 

2. Comparison of the 880-nm data to contemporaneous dust opti- 
cal depth measurements from Mini-TES shows that changes in 
the Mini-TES/880-nm ratio frequently occur with the onset 
and decay of dust storms. The values are consistent with varia- 
tions in the mean radius of dust aerosols from near 1 pm during 
northern summer to over 2 pm at the onset of a storm, and 
being generally near the reported values of 1.4 pm (from obser- 
vations near the beginning of the mission). 
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